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A mathematical model of the process of blowing high-temperature products of fuel combustion through mol¬ 
ten glass in the melting tank of a glassmaking furnace via side tuyeres placed below melt level is presented. 
Blowing is conducted in two regimes — bubble and jet. 
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Gas can be blown through molten glass via a tuyeres, lo¬ 
cated below the melt level, in two regimes — bubble and jet. 
In the first case the gas at the tip of a tuyere forms a gas ca¬ 
vity which is held on the tip by surface tension. The 
expulsive buoyancy force strives to detach the cavity. As the 
volume of the gas cavity increases the buoyancy force in¬ 
creases, while the surface tension remains unchanged. In the 
limiting case these forces are equal to one another. As it con¬ 
tinues to grow the cavity detaches from the tip of the tuyere, 
and in the process a gas bubble forms and floats upward in 
the melt tank. 

In the second case, where a gas jet flows into the molten 
glass, the gas at the exit from the tuyere forms a stable gas 
cavity, the flow in which forms a gas jet. Individual bubbles 
detach from the surface of the gas jet, while the cavity itself 
is continually sits on the tip of the tuyere. 

In the present work we examine the case where 
high-temperature products of fuel combustion are blown into 
the molten glass in the melt tank through a side tuyere placed 
below the melt level. 

The following assumptions are made to construct a 
model of the blowing process: 

- the gas and liquid are chemically neutral, isothermal 
and incompressible; 

- the molten glass is a homogeneous viscous liquid; 

- before detaching from the tip of the tuyere a gas bub¬ 
ble is symmetric relative the longitudinal axis of the tuyere; 

- the gas jet is symmetric relative to the longitudinal 
axis of the tuyere over the entire distance to breakup into 
bubbles; 
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- the minimum static gas pressure at the exit from 
tuyere equals the static pressure in the furnace space above 
the molten glass. 

The assumption of the chemical neutrality of the gas and 
liquid is adopted because in a physical model it is impossible 
to simulate simultaneously the hydrodynamics and mass 
transfer accompanied by a chemical reaction. 

The gas and liquid are assumed to be incompressible be¬ 
cause the pressure in the bubbling layer varies within negligi¬ 
ble limits while the velocity of these phases is relatively low. 

For a number of technological apparatuses with a bub¬ 
bling layer the temperatures of the gases flowing into and out 
of the layer differ by not more than 10 - 15%, i.e., the gas 
and liquid in the layer are practically isothermal. In cases 
where a cold gas is blown into a tank of molten glass the 
temperature of the gas becomes equal to the temperature of 
the molten glass within less than two tuyere lengths. In other 
words, even in this case the gas and liquid in most of the inte¬ 
rior volume of the layer interact as isothermal media. On the 
initial section, where the gas temperature increases, the ve¬ 
locity of the gas increases, which should intensify mixing of 
the liquid in a zone near the tuyere. In this connection a study 
of the efficiency of heat and mass transfer using an isother¬ 
mal hydraulic model corresponds to the least favorable con¬ 
ditions for the process to proceed. 

It is known that the mass content of the disperse phase 
(solid particles or drops of insoluble liquid) in a bubbling 
layer does not exceed 10% 3 [1]. Impurity content to 20% of 
the total mass of the melt has no appreciable effect on the hy¬ 
drodynamics of the bubbling layer [2], Therefore, it can be 
assumed that the viscous liquid (molten glass) is homoge¬ 
neous in the melt tank. 


3 Here and below the content by weight. 
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When gas is blown through a side tuyere in the bubble i.e., when 
regime a growing gas bubble moves upward relative to the 
longitudinal axis of the tuyere. At the moment of detachment 
the bubble is located above the top edge of the tuyere. As a 
result the gas pressure inside the bubble drops by the amount 
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where g is the acceleration of free fall, m/sec 2 ; p gas is the 
density of the gas, g/m 3 ; and, d 0 is the diameter of the open¬ 
ing of the tuyere, m. 

The gas pressure inside the bubble is given by 
Pg ftPliq Yq J 


Taking p liq = 2300 kg/m 3 , A liq =1.0m and p gas = 
1.29 kg/m 3 we find that this assumption definitely holds for 
(w gas ) max >189 m/sec. Since the velocity obtained is close to 
the average outflow velocity of the gas in most furnaces with 
a bubbling layer, this assumption is correct. 

We now represent the mathematical model of the hydro¬ 
dynamics of the bubbling layer in the following form: 


where p liq is the density of the liquid (molten glass), kg/m 3 , 
and /i| iq is the height of the layer of liquid above the axis of 
the tuyere, m. 

Since p g is at least three orders of magnitude greater than 
A p, any change of the gas pressure inside a bubble can be ne¬ 
glected and it can be assumed that as it grows a bubble re¬ 
mains symmetric relative to the longitudinal axis of the 
tuyere. 

The sum of the dynamic and static gas pressures at the 
exit from the tuyere is a constant, equal to the sum of the 
counter pressure of the molten glass and turbulent pulsations. 
In accordance with [1], the energy of the turbulent pulsations 
is practically a constant; during blowing the ratio between 
the dynamic and state pressures varies correspondingly. At 
the moment the gas cavity detaches from the tip of tuyere the 
liquid covers the entrance opening of the tuyere, and the gas 
stops for an instant. At this moment the static pressure of the 
gas has its maximum value equal to the stopping pressure: 

(pTr ax =Po + SP liq h liq , 

where p 0 is the static pressure on the melt surface, Pa, and 
h u is the height of the fined (gas-free) liquid above the axis 
of the tuyere, m. 

At the next instant a gas cavity starts to form, under the 
pressure of the gas, in front of the tuyere tip. The gas begins 
to move, and as the cavity grows in size its velocity in¬ 
creases. The total pressure developed by the gas at the tuyere 
exit is 
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Before the next bubble detaches from the tip of the 
tuyere the gas velocity w gas assumes its maximum value and 
the static gas pressure its minimum value. 

The assumption that the minimum static gas pressure at 
the tuyere exit equals the static gas pressure in the furnace 
space above the molten glass (p| as ) min ~ p 0 holds if 
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where w gas , w liq are, respectively, the velocity of the gas and 
liquid, m/sec; k is the mass density of the external mass 
force (gravity), N/kg; p gas , p liq are the density of the gas and 
liquid, respectively, kg/m 3 ; p |iq is the pressure of the liquid at 
the surface of a bubble. Pa; p ga& is the gas pressure inside the 
bubble, Pa; grad p is the gradient of the pressure, Pa/m; v„ as , 
v u are, respectively, the kinematic coefficients of viscosity 

q q2 q2 q2 

of the gas and liquid, m 2 /sec; V 2 =—- + is the 

dx 2 dy 1 dz 2 

. , du dv dw . , 

Laplace operator; div(w) = — + — + — is the divergence 

dx dy dz 

of the velocity vector; w liq _ gas is the relative velocity of lav¬ 
ing of the gas bubble by the molten glass, m/sec; g is the 
acceleration of free fall, m/sec 2 ; /? gas , /A are the pressure of 

the gas and liquid on the phase interface. Pa; p gas , p liq are the 
dynamic coefficient of viscosity of the gas and liquid, re¬ 
spectively, Pa • sec; o is the interphase tension at the gas bub¬ 
ble - liquid interface, N/m; R v R 2 are the principal radii of 
curvature of the surface of the rising bubble (the bubble is as¬ 
sumed to be an ellipsoid), m; and, w h is the rise velocity of 
the bubble, m/sec. 

The uniqueness conditions are formulated separately for 
the bubble and jet regimes of outflow: 
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for the bubble regime: 
- for t > t n 


- for X=0 

- for Z = 0 
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for the jet regime: 
- for t = tn 
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where t Q is the moment of onset of bubble growth or a transi¬ 
tion of the bubble into the jet outflow regime; d b is the bub¬ 
ble diameter, m; d-. is the diameter of the transverse section 
of the jet, m; iv 0 is the average flow velocity of the gas in the 
exit section of the tuyere, m/sec; w cr is the velocity at which 
the jet outflow regime commences, m/sec; /j et is the length of 
the jet in the melt, m; c, is the displacement along the Z axis 
of the fluctuating gas jet - molten glass interface; X is the 
wavelength of these oscillations, m; b is the amplitude value 
of the fluctuating acceleration, acting along the vertical axis 
Z and due to the oscillations of the jet - molten glass inter¬ 
face, m/sec 2 ; and, eo is the circular frequency of the oscilla¬ 
tions, rad/ sec. 

We shall assume the outflow to be quasistationary. Then, 
under the uniqueness conditions the following are also pre¬ 
scribed: 

- the geometric parameters of the furnace (interaxial 
distance of the tuyeres / tu , m; diameter of the outflow open¬ 
ing of a tuyere d 0 , m; height of the layer of fined molten glass 
above the axis of the tuyere K q , m); 

- the physical parameters of the molten glass and gas 
(p, kg/m 3 ; p, Pa • sec; a, N/m); 

- the gas velocity at the exit from a tuyere w 0 , m/sec; 


- the sound velocity in the gas a, m/sec. 

The equations (1) - (4) of the system (I) are the equa¬ 
tions of motion and continuity of the gas and liquid. The 
equations (5), (6) are the equations of motion and continuity 
of the liquid flowing around a gas bubble, conventionally as¬ 
sumed to be stationary. The equation (7) expresses the condi¬ 
tion that the normal stresses on the surface of a bubble taking 
account of the surface tension and neglecting the viscosity of 
the gas are equal. The equation (8) expresses the absence of 
phase slip at the surface of a bubble. 

The physical meaning of the system (II) is as follows. In 
the bubble regime of outflow, at the moment a bubble de¬ 
taches liquid momentarily covers the opening, in the exit sec¬ 
tion of a tuyere, the gas moving along the delivery channel is 
suddenly decelerated right up to complete stoppage, and its 
kinetic energy is converted into a static head. This assumes 
that /; liq » d 0 . The symbol X denotes the distance along the 
longitudinal axis of the tuyere at which gas decelerates from 
the velocity w 0 to zero and t is the time during which this de¬ 
celeration occurs. 

The system (III) describes the mechanism of the suppres¬ 
sion of the expulsive force on the surface of the gas jet due to 
the turbulent pulsations of the gas. 

The condition that the dynamic pressure of the gas flow 
equals the static pressure of the molten glass in the melt tank 
of the furnace always holds in the jet outflow regime at the 
end of the jet. It is expressed by the last equation of the sys¬ 
tem (III). 

The mathematical model presented above is the theoreti¬ 
cal basis for developing the methods and apparatus for physi¬ 
cal modeling of mixing processes in molten glass through 
which a gas jet, located on the side wall of the furnace below 
the melt level, is blown. 
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